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FOREWORD 


This  report  was  prepared  in-house  under  Project  No.  1347,  "Structural  Testing  of 
Flight  Vehicles,”  Task  No.  134703,  "Structural  Testing  Criteria.”  The  work  was 
conducted  by  the  Experimental  Mechanics  Branch,  Structures  Division,  Air  Force  Flight 
Dynamics  Laboratory,  with  Mr.  W.  R,  Johnston  acting  as  task  engineer. 
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ABSTRACT 


This  study  was  conducted  to  determine  how  accurately  aerodynamic  heating  is 
simulated  in  the  laboratory  with  radiant  heat  lamps  under  computer  control.  A  study 
of  aerodynamic  heating  theory  resulted  in  the  preparation  of  a  digital  computer  program 
for  calculating  the  parameters  required  as  input  data  for  the  heat  rate  computer.  The 
heat  rate  computer  is  used  to  control  the  heat  input  to  an  actual  flight  test  specimen. 
Laboratory  test  results  are  compared  with  flight  measured  data  to  determine  the  degree 
of  simulation.  Test  results  indicate  that  aerodynamic  heating  may  be  adequately 
simulated  up  to  600°F.  Laboratory  simulation  to  higher  temperatures  are  being  studied 
for  future  programs. 


This  technical  documentary  report  has  been  reviewed  and  is  approved. 


HOLLAND  B.  LOWNDES,  JR. 
Actg  Chief,  Structures  Division 
AF  Flight  Dynamics  Laboratory 
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SYMBOLS 

SYMBOL 

DEFINITION 

UNITS 

A 

Area 

Fta 

B 

Emissivity  factor  (ecr) 

BTU/hr-ft3  -°R' 

B-h 

B-h  SCALE  FACTOR  potentiometer  setting 

— 

c 

Specific  heat  of  the  skin  material 

BTU/lb-°R 

c 

p 

Specific  heat  of  air 

BTU/Ib-°R 

c 

put 

C  SCALE  FACTOR  potentiometer  setting 

— 

^switch 

C  SCALE  FACTOR  switch  position 

— 

E 

Voltage 

Volts 

Ae 

Error  signal 

Volts 

El  gain 

Inherent  gain  associated  with  the  El  RANGE 
switch  position 

— 

G 

Gain 

— 

G* 

Nocturnal  irradiation 

BTU/hr-ft3 

G 

s 

Solar  irradiation 

BTU/hr-ft3 

g 

Acceleration  due  to  gravity 

Ft/sec3 

h 

Convective  heat  transfer  coefficient 

BTU/hr-ft3 -°R 

H 

aw 

Adiabatic  wall  enthalpy 

BTU/lb 

H 

s 

Skin  enthalpy 

BTU/lb 

I 

Current 

Amperes 

J 

Conversion  factor  (778) 

Ft-lb/BTU 

kd 

Demand  gain 

— 

K 

pot 

K  SCALE  FACTOR  potentiometer  setting 

— 

v 
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SYMBOL 
P.B.  Offset 

% 

^loss 


r 

T 


aw 


'  stag 


sp 


st 


s  range 


t 

V 
w 

Y 


a 


SYMBOLS  (Continued) 
DEFINITION 
Proportional  band  offset 
Heat  input  from  equipment 

Total  heat  losses  due  to  internal 
convection,  conduction,  and  radiation 

Recovery  factor 

Adiabatic  wall  temperature 

Local  free  stream  temperature 

Effective  temperature  of  space 

Skin  temperature 

Stagnation  temperature 

Programmed  skin  temperature 

Test  skin  temperature 

T  RANGE  switch  position 
s 

Time 

Velocity 

Specific  weight  of  the  skin  material 
Conversion  factor 
Nocturnal  absorbtivity 

Solar  absorbtivity 

Skin  emissivity 


UNITS 

°R 

BTU/hr-ft3 

BTU/hr-ft3 


°R 

°R 

°R 

°R 

°R 

°R 

°R 

°R 

Hr 

Ft/sec 

Lb/ft3 

BTU/hr-watt 


vi 


FDLTDR  64-112 
Part  I 


SYMBOLS  (Continued) 

SYMBOL  DEFINITION 

n  Efficiency 

a  Stefan -Boltzmann  constant  (1.73  x  10-9) 

r  Skin  thickness 


UNITS 


BTU/hr-fta-°R'* 

Ft 
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COMPUTER  PROGRAM  SYMBOLS 

SYMBOL 

DEFINITION 

UNITS 

H 

Altitude 

Feet 

VA 

Velocity  of  the  vehicle 

Ft/sec 

THETA 

Time 

Sec 

TS 

Skin  temperature 

Or 

DTS 

Derivative  of  skin  temperature  with  respect 
to  time,  dT/d0 

°R/sec 

XMN 

Free  stream  Mach  number 

— 

T2 

Temperature  of  the  air  aft  of  the  shock  wave 

°R 

RN2 

Reynolds  number  aft  of  the  shock  wave 

— 

TTS 

Temperature  values  for  table  used  to  find 
skin  emissivity  and  specific  heat 

“R 

ES 

Emissivity  of  the  skin  material 

— 

CPS 

Specific  heat  of  the  skin  material 

BTU/lb-°R 

DTSI 

Derivative  of  skin  temperature  with  respect 
to  the  time  used  in  the  Runge-Kutta  calculations 

°R/sec 

ALFD 

Angle  of  attack 

Degrees 

TAW 

Adiabatic  wall  temperature 

“R 

HTC 

Aerodynamic  heat  transfer  coefficient 

BTU/hr-ft3  -°R 

Q 

Total  heat  transferred 

BTU/hr-ft3 

TAWP 

Adiabatic  wall  temperature  scaled  for  punching 
on  tape  for  use  on  the  heat  rate  computer 

— 

HP 

Aerodynamic  heat  transfer  coefficient  scaled 
for  punching  on  tape  for  use  on  the  heat  rate 
computer 

TSP 

Skin  temperature  scaled  for  punching  on  tape 
for  use  on  the  heat  rate  computer 

— 

K 

Number  of  points 

»  -  « 
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COMPUTER  PROGRAM  SYMBOLS  (Continued) 


SYMBOL 

DEFINITION 

UNITS 

PRP 

Prandtl  Number  raised  to  the  -2/3  power 

— 

RNCR 

Critical  Reynolds  number 

— 

D7 

Wedge  angle 

Degrees 

ELH 

Distance  from  leading  edge  to  the  point  in 
question 

Feet 

RHOS 

Density  of  the  skin  material 

Lbs /ft3 

DELTAS 

Skin  thickness 

Feet 

DELTAH 

Time  interval  used  in  Runge-Kutta  calculations 

Sec 

ESA 

Average  material  emissivity 

— 

HGP 

Altitude 

Meters 

Al,  A2,  A3 

Constants  for  atmospheric  calculations 

— 

TB 

Atmospheric  temperature  constant 

°R 

PB 

Atmospheric  pressure  constant 

Lbs /ft3 

RHOB 

Atmospheric  density  constant 

Slugs/ft3 

HB 

Altitude  constant 

Meters 

P 

Local  atmospheric  pressure 

Lbs/fta 

RHO 

Local  atmospheric  density 

Slugs/ft3 

TM 

Local  molecular-scale  temperature  of  the 
atmosphere 

°R 

VS 

Local  speed  of  sound 

Ft /sec 

VIS 

Local  kinematic  viscosity  of  the  atmosphere 

Ft3  /sec 

T 

Local  temperature  of  the  atmosphere 

°R 

A,  B,  C,  D 

Constants  for  finding  the  atmospheric 
temperature 

- - 

ALFHD 

Wedge  angle  plus  angle  of  attack  (D7+ALFD) 

Degrees 
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COMPUTER  PROGRAM  SYMBOLS  (Continued) 


SYMBOL 

DEFINITION 

UNITS 

P2 

Pressure  aft  of  the  shock  wave 

Lbs/ft3 

RH02 

Density  of  the  air  aft  of  the  shock  wave 

Slugs/ft3 

V2 

Velocity  aft  of  the  shock  wave 

Ft/sec 

ALFHR 

(D7  +  ALFDJ/57.3 

Radians 

CP 

Pressure  coefficient 

— 

XMNN 

Mach  number  normal  to  the  shock  wave 

— 

BETA 

Shock  wave  angle 

Radians 

XMN2 

Mach  number  aft  of  the  shock  wave 

— 

BMU2 

Coefficient  of  viscosity  aft  of  the  shock  wave 

Lb/ft-sec 

CH,  YH,  RH 

Constants  depending  on  Reynolds  Number  (RH  is 
the  recovery  factor) 

— 

Dl,  D2,  D3 

Constants  used  in  calculating  enthalpies 

— 

H2 

Enthalpy  aft  of  the  shock  wave 

BTU/lb 

HAW 

Adiabatic  wall  enthalpy 

BTU/lb 

HS 

Skin  enthalpy 

BTU/lb 

HR 

Reference  enthalpy 

BTU/lb 

D4,  D5,  D6 

Constants  used  in  calculating  temperatures 

— 

TH 

Temperature  based  on  reference  enthalpy 

°R 

BMUH 

Coefficient  of  viscosity  based  on  reference  enthalpy 

Lb/ft-sec 

RHOH 

Density  based  on  reference  enthalpy 

Lb/ft3 

RNH 

Reynolds  Number  based  on  reference  enthalpy 

— 

SBC 

Stefan-Boltzmann  constant  (4.758  x  10~13) 

BTU/ sec-ft3  -  °R 

TSA 

Temperature  used  in  Runge-Kutta  calculations 

°R 
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COMPUTER  PROGRAM  SYMBOLS  (Continued) 


SYMBOL 

DEFINITION 

UNITS 

DTA,  DTB,  DTC 

Derivative  of  temperature  with  respect  to  time 
used  in  Runge-Kutta  calculations 

°R/sec 

CPA 

Specific  heat  of  air 

BTU/lb-°R 

BC 

Average  emissivity  multiplied  by  the  Stefan- 
Boltzmann  constant 

BTU/sec-fta-°R4 

HMAX 

Maximum  aerodynamic  heat  transfer  coefficient 

BTU/hr-ft®  -°R 

TAWMAX 

Maximum  adiabatic  wall  temperature 

°R 

BH 

Calculated  setting  for  the  heat  rate  computer 

B-h  potentiometer 

— 

DK 

Demand  gain  obtained  in  scaling  for  the  heat 
rate  computer 

— 

TAWR 

Adiabatic  wall  temperature  range  for  the  heat 
rate  computer 

°R 

TSMAX 

Maximum  skin  temperature 

°R 

TSR 

Skin  temperature  range  for  the  heat  rate 
computer 

°R 
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INTRODUCTION 


This  program  was  conducted  to  determine  if  elevated  temperature  test  techniques 
used  in  the  Structures  Test  Facility  of  the  Air  Force  Flight  Dynamics  Laboratory 
provide  an  acceptable  simulation  of  aerodynamic  heating.  A  background  study  of 
aerodynamic  heating  theory  was  performed  and  a  digital  computer  program  was 
written  to  provide  the  necessary  information  for  use  on  a  heat  control  computer. 

Thermal  simulation  tests  were  conducted  on  a  flight  test  article,  and  the  laboratory 
test  temperatures  were  compared  with  the  flight  test  temperatures  to  determine  the 
degree  of  simulation. 

AERODYNAMIC  HEATING 

When  a  flight  vehicle  moves  through  the  earth’s  atmosphere,  the  relative  motion 
between  the  vehicle  and  the  air  creates  a  flow  boundary  layer  (Reference  I).  The 
frictional  forces  between  the  vehicle  and  the  boundary  layer  convert  the  energy  of 
the  moving  air  into  heat.  The  amount  of  heat  transferred  to  the  vehicle  surface  is 
dependent  upon  the  difference  between  the  boundary  layer  temperature  and  the  skin 
temperature,  and  the  coefficient  of  heat  transfer.  The  rate  of  change  of  the  skin 
temperature  depends  on  the  rate  of  heat  exchange  and  the  heat  capacity  of  the  skin 
structure.  If  steady -state  aerodynamic  conditions  are  maintained  for  a  sufficient  time, 
a  condition  will  result  wherein  aerodynamic  heat  input  equals  the  total  heat  losses,  thus 
resulting  in  a  constant  surface  temperature  . 

Some  of  the  heat  generated  by  the  aerodynamic  conditions  is  dissipated,  or  lost,  from 
the  vehicle’s  surface  by  radiation.  At  low  velocities  the  skin  temperatures  are  low  and 
the  resultant  radiation  losses  are  an  insignificant  portion  of  the  total  heat  transferred. 

At  large  Mach  numbers,  and  therefore  high  temperatures,  the  amount  of  heat  radiated 
by  the  surface  may  be  a  large  percentage  of  the  heat  transferred. 

The  aerodynamic  heating  properties  are  based  on  the  local  free  stream  conditions 
at  the  point  of  the  vehicle  under  consideration.  At  velocities  greater  than  the  speed  of 
sound,  shock  and  expansion  areas  are  produced  that  create  local  conditions  different 
from  the  surrounding  ambient  atmospheric  conditions.  The  changes  that  take  place  are 
dependent  upon  the  shape  of  the  vehicle.  All  of  the  calculations  in  this  study  are  based  on 
a  wedge-shaped  body.  Although  a  large  amqunt  of  the  information  in  this  report  is 
applicable  to  any  shape,  many  of  the  calculations  in  the  digital  computer  program 
would  require  changing  for  use  with  any  shape  other  than  a  wedge. 


Manuscript  released  by  author  7  July  1964  for  publication  as  an  RTD  Technical 
Documentary  Report. 
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THE  HEAT  BALANCE  EQUATION 


In  order  to  investigate  aerodynamic  heating  thououghly,  all  possible  heat  inputs  and 
losses  must  be  recognized.  For  illustrative  purposes,  consider  a  flat  plate  insulated 
on  all  surfaces  except  the  heated  surface  and  subjected  to  an  aerodynamic  heating 
environment.  The  immediate  environmental  temperature  will  depend  only  on  the  vehicle 
altitude  and  velocity.  The  stagnation  temperature  is  expressed  by  the  following  equation: 


The  temperature  of  the  air  stream  next  to  the  surface  (the  adiabatic  wall  temperature) 
is  related  to  the  ambient  temperature  and  the  stagnation  temperature  by  the  following 
expression: 


T 


2g  JCp 


The  amount  of  heat  that  reaches  the  surface  of  the  vehicle  is  dependent  upon  the  value 
of  the  convective  heat  transfer  coefficient.  This  heat  is  represented  by  the  expression 
h(T  -  Tg).  Additional  heat  may  be  obtained  by  irradiation  from  surrounding  bodies 

such  as  the  earth  and  the  sun.  If  the  vehicle  is  radiating  to  space,  the  surface  will  lose 
heat  according  to  the  term  ctrT  The  heat  capacity  of  the  plate  per  unit  time  can  be 

d  "|  ® 

expressed  as  cwr  —  .  The  previous  terms  can  be  combined  to  form  the  following  the 
balance  equation: 


-  TS  >  +  a@  g{3 


It  must  be  emphasized  that  this  expression  represents  the  idealized  case  of  a  uniformly 
heated  insulated  plate  with  a  unit  area,  constant  thickness,  and  uniform  material 
properties. 

Since  the  idealized  case  seldom  exists  in  actual  flight,  the  basic  heat  balance  equation 
must  be  modified  to  include  other  sources  of  heat  inputs  and  losses.  Additional 
heat  inputs  may  be  obtained  from  engines,  electronic  devices,  electric  motors,  and 
other  heat  generating  equipment.  Heat  may  be  lost  by  conduction  to  adjacent  structure, 
internal  convection  and  radiation,  and  through  active  cooling  equipment.  The  heat 
balance  equation  can  now  be  written  as  follows: 


h{To*-Ts  1  +  af3  Gf3  +  as  Gs  +  -  ecr  TS  -^loss 


where  q  represents  additional  heat  inputs,  and  q  represents  additional  heat  losses. 
0  lOSS 
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During  accelerated  or  decelerated  atmospheric  flight  some  of  the  terms  become 
relatively  insignificant.  For  instance,  solar  and  nocturnal  irradiation  is  significant  only 
in  space  operations.  Similarly,  the  heat  generated  by  equipment  is  often  a  very  small 
percentage  of  the  heat  generated  by  aerodynamic  heating.  By  eliminating  these  insignifi¬ 
cant  factors  the  heat  balance  equation  may  be  expressed  as  follows: 


This  equation  is  the  basis  for  the  closed  loop,  computer  controlled  aerodynamic  heating 
simulation  tests. 


THE  TEST  PROGRAM 

When  conducting  elevated  temperature  tests,  the  philosophy  is  to  simulate  the  external 
thermal  environment  and  allow  the  structure  to  distribute  the  heat  as  it  would  in  actual 
flight.  This  method  is  in  contrast  to  the  philosophy  of  adjusting  the  external  test 
environment  until  some  predicted  internal  temperature  distribution  is  obtained.  In 
addition,  the  specimen  is  used  as  a  calorimeter,  and  the  surface  temperature  is  fed 
back  from  each  control  area  to  the  computer  which  calculates  any  power  change 
required.  This  instrumentation  and  control  feature  contrasts  to  the  method  in  which 
heat  flux  gages  were  used  by  some  French  and  British  investigators. 

It  is  obvious  that  radiant  heating  is  different  from  convective  heating  and  does  not 
duplicate  the  mass  flow  encountered  by  a  vehicle  moving  through  the  earth’s  atmosphere. 
The  objectives  of  this  program  were  to  determine  the  validity  of  the  laboratory’s  test 
philosophy  and  the  accuracy  of  aerodynamic  heating  simulation  using  radiant  heat 
lamps  under  computer  control.  If  the  temperature  distribution  obtained  with  radiant 
heating  is  the  same  as  that  obtained  during  flight,  then  it  may  be  assumed  that  the 
aerodynamic  thermal  environment  has  been  adequately  simulated.  In  this  study  the 
laboratory  test  temperatures  were  compared  with  flight  measured  temperatures  to 
determine  the  accuracy  of  laboratory  simulation. 

^  This  study  is  also  concerned  with  the  adequacy  of  the  design  equations  for  the  General 
Electric  Heat  Rate  Computer,  which  is  used  for  control  of  most  elevated  temperature 
tests  conducted  in  this  laboratory.  Of  particular  interest  is  the  computer  limitation  that 
constants  be  used  for  temperature  dependent  variables.  For  instance,  the  specimen 
emissivity  and  specific  heat  are  known  to  vary  with  temperature,  but  a  constant  value 
must  be  used  in  setting  up  the  computer  (Reference  2).  Similarly,  the  efficiency  of  a 
test  setup  will  vary  as  the  temperature  and  temperature  rise  rate  change,  but  again 
the  computer  requires  that  the  efficiency  be  considered  a  constant.  In  addition,  since  the 
computer’s  design  equations  were  based  on  an  idealized  insulated  plate,  information 
was  desired  on  what  effect  adjacent  heat  sources,  heat  sinks,  and  various  structural 
members  would  have  on  the  computer's  control  characteristics.  The  effects  these  factors 
have  on  thermal  simulation  tests  are  important  for  determining  design  requirements  for 
future  heat  control  computers. 
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Two  test  specimens,  with  complete  flight  data,  were  available.  A  nose  cone  specimen 
(a  polished  stainless  steel  truncated  cone  shell)  was  fabricated  to  be  a  duplicate  of  a 
Viking  10  rocket  nose  cone.  Information  for  nose  cone  construction  and  test  correlation 
was  obtained  from  NRL  Report  4531.  "Flight  Measurements  of  Aerodynamic  Heating  and 
Boundary -Layer  Transition  on  the  Viking  10  Nose  Cone”  (Reference  3).  Care  was  taken 
in  the  laboratory  setup,  shown  in  Figures  1  and  2,  to  duplicate  the  flight  conditions  as 
closely  as  possible.  Twenty-five  thermocouples  were  installed  on  the  specimen, 
nineteen  internal  data  thermocouples  and  six  external  control  thermocouples.  The  nose 
cone  was  divided  into  three  thermal  control  areas,  with  maximum  temperatures  of 
approximately  300°F.  Because  of  low  temperatures  and  resultant  low  power  require¬ 
ments,  the  nose  cone  was  abandoned  in  favor  of  the  X-7A  wedge-shaped  wing. 

The  X-7A  wing,  shown  in  Figure  3,  was  originally  designed  to  investigate  the  effects 
of  transient  aerodynamic  heating  on  aircraft  structures.  The  leading  edge  was  ground  to 
produce  a  knife-edge  wedge  with  an  included  angle  of  nine  degrees,  insulative  materials 
were  applied  to  some  areas,  boundary  layer  trips  were  installed  to  produce  turbulent 
flow,  and  some  areas  were  etched  to  provide  varying  degrees  of  surface  roughness. 

In  addition  to  the  skin  thermocouples,  thermocouples  were  installed  on  structural 
members  to  determine  the  internal  temperature  distribution.  Construction  details  and 
flight  test  data  are  available  in  Reference  4. 

The  X-7A  wing  test  setup  is  shown  in  Figure  4.  The  wing  was  placed  between  two 
polished  aluminum  reflectors,  on  which  were  installed  480  T-3  heating  lamps  (1000- 
watt).  The  wing  was  divided  into  eight  thermal  control  areas,  four  for  the  top  surface 
and  four  for  the  lower  surface.  The  maximum  skin  temperature  required  was  approxi¬ 
mately  600°F.  Temperatures  were  measured  using  the  actual  flight  test  thermocouples 
which  consisted  of  41  skin  and  6  structural  thermocouples  installed  inside  the  wing. 

Eight  additional  thermocouples  were  installed  on  the  wing  external  surfaces  for  computer 
control  purposes.  The  temperature  was  recorded,  reduced  and  displayed  "on-line”  by 
means  of  ETSTF  HSDAPS  (Elevated  Temperature  Structural  Test  Facility's  High 
Speed  Data  Acquisition  and  Processing  System). 

CONCLUSIONS 

Satisfactory  results  were  obtained  in  all  tests  controlled  by  the  General  Electric  Heat 
Rate  Computer.  Laboratory  temperatures  were  within  5  percent  of  the  flight  recorded 
temperatures.  This  is  an  acceptable  accuracy  when  it  is  considered  that  the  information 
originated  with  only  the  flight  profile  data,  the  calculations  were  performed  by  a  digital 
computer,  and  the  final  test  was  controlled  by  a  closed-loop  computer. 

The  program  test  results  have  shown  that  aerudynamic  heating  can  be  accurately 
simulated  up  to  600°F  with  radiant  heat  lamps  under  computer  control.  However,  some 
of  the  areas  of  original  interest  still  require  more  study.  For  example,  the  temperatures 
reached  to  date  are  too  low  to  reveal  any  significant  information  concerning  the  effect 
of  using  constants  in  place  of  temperature  dependent  variables.  Also,  in  the  X-7A  wing 
test,  no  attempt  was  made  to  cool  the  wing  to  the  sub-zero  temperatures  experienced 
at  the  beginning  of  the  test  flight.  More  accurate  results  would  probably  have  been 
obtained  if  the  wing  had  been  initially  cooled. 
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Figure  1*  Front  View  of  Nose  Cone  Test  Setup 
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Figure  Zt  Rear  View  of  Nose  Cone  Test  Setup 
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Figure  3,  X-7A  Wedge-Shaped  Wing 
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Figure  4,  X-?A  Wing  Test  Setup 


FDL  TDR  64-112 
Part  I 

Future  programs  are  planned  to  study  simulation  techniques  to  higher  temperatures 
where  radiation  losses  and  temperature  dependent  variables  become  more  significant. 
The  test  programs  will  be  controlled  by  a  more  versatile  heat  control  computer  and 
use  a  more  elaborate  test  setup  to  duplicate  the  total  flight  environment  as  much  as 
possible.  It  is  anticipated  that  the  information  obtained  will  enable  this  laboratory  to 
confidently  and  accurately  simulate  aerodynamic  heating  on  any  aerospace  vehicle. 
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APPENDIX  I 

IBM  DIGITAL  COMPUTER  PROGRAM 

More  efficient  performance  of  the  many  calculations  required  in  this  study  were 
achieved  by  the  digital  computer  program  written  for  use  on  the  IBM  7094  computer. 
The  program  is  in  the  engineers  computer  language,  FORTRAN,  and  is  based  on  the 
aerodynamic  heating  subprogram  given  in  Reference  5.  The  program  shown  as 
Table  I,  has  been  designed  to  provide  information  peculiar  to  the  operation  of  the 
General  Electric  Heat  Rate  Computer  used  in  this  laboratory,  but  it  still  retains  its 
FORTRAN  style  and  may  be  easily  utilized  by  anyone  with  a  knowledge  of  FORTRAN 
programming. 

The  calculations  are  based  on  the  aerodynamic  heating  theory  for  a  thin-skinned 
wedge  at  an  angle  of  attack.  By  ingoring  conduction  into  the  structure,  the  heat  energy 
stored  in  the  skin  is  considered  as  the  difference  between  the  aerodynamic  heat  input 
and  the  heat  radiated  to  space.  From  this  theory  the  following  differential  equation  is 
derived: 


111 

dt 


c w  TC„ 


€<x 


TT 'T.‘  -V  > 


where  c  and  e  are  properties  of  the  skin  material  and  surface  condition,  and  h,  C  ,  and 

P 

H  are  properties  of  the  air  flowing  over  the  point  under  consideration.  All  of  these 

5 

items  are  functions  of  the  skin  temperature.  Since  the  resulting  differential  equation  has 
non-linear  coefficients,  the  Runge-Kutta  numerical  method  is  used  to  obtain  the  solution. 
The  calculations  for  determining  the  heat  transfer  properties  are  based  upon  a  reference 
enthalpy  method,  and  calculations  for  determining  atmospheric  conditions  and  local  flow 
conditions  are  based  upon  the  1959  ARDC  Model  Atmosphere. 

As  it  is  written  the  digital  computer  program  incorporates  the  following  features: 

1.  Up  to  300  time  points  for  each  location  on  the  wedge 

2.  Altitudes  up  to  700,000  meters 

3.  Atmospheric  properties  calculated  according  to  the  1959  ARDC  Model 
Atmosphere 

4.  Provisions  for  the  reading  of  multiple  sets  of  data  with  easy  variation  of 
Prandtl  Number,  critical  Reynolds  Number,  material,  and  location  of  each 
point  in  question 

5.  Provisions  for  a  linear  interpolation  table  for  variation  of  material  specific 
heat  and  emissivity  with  skin  temperature 

6.  Runge-Kutta  method  of  solving  the  non-linear  differential  equation 
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TABLE  I 

IBM  DIGITAL  COMPUTER  PROGRAM 

C  SKIN  TEMPERATURE  PROGRAM 

DIMENSION  H ( 300  )  *  V  A [  3  00  )  »  T  HET  A (  300  )  *  T  ^  0 1  )  »PT6(30n),xMN[3nO)  iTJIi 

100 ) *RN2( 300) ,TTS( 20  I ,ES( 20) ,CPS( 20  I .OTSU  4) *ALFD(300 I »TAWt  300  I »HTC 
2  t  300  )  >0(300) >TAWP I  300) *HP( 300 ) >TSP ( 300  > 

A  READINPUTTAPE?,1  ,K 

1  FORMAT  (13) 

RF AO  INPUT  TAPE  2 > 2 > ( THF T A < M ) > M  =  1 > K  ) 

RFAD  INPUT  TAPE  2,2,  (H(M)  »  M  =  i  ,  K  ] 

READ  INPUT  TAPE  2 . 2 > ( V A ( M ) > M.=  1 » K  ) 

RF AO  INPUT  TAPE  2,2t( ALED(M) ,M=1 »K) 

REAP  INPUT  TAPE  2 >2 , ( TT S ( M I , M= 1 , 20 ) 

READ  INPUT  TAPE  2 > 2 , I  CPS ( M ) > m= 1 > 2 0 ) 

READ  INPUT  TAPE  2 , 2  ,  ( ES ( M > » M  =  1 » 20 }  . 

2  FORMAT  ( 5 El  3 , 4 ) 

5  RFADINPtJTTAPE?»3,T5(  1  )  ,  PRP  >  RNCR  *  07,  ELH»RHOS>  DELTAS  »DELTAH*ESA 

3  FORMAT  t  E  1 1  • 4 ) 

DO  200  N  =  1  *  K 

HGP  =  .3048*H[N ) / t 1 ,  +  . 3048*H(N) /63S6766.  ) 


IF 

( HGP-7  00000 . I 

10,20,213 

10 

IF 

(HGP-200000. ) 

11,30,20 

11 

IF 

( HGP- 1 70000. ) 

1 2  ,40,30 

12 

IF 

(HGP-160000. ) 

13,50,40 

1  ? 

IF 

(HGP-1 06000. ) 

14,60,50 

I  A 

IF 

(  HGP-9 0000 . ) 

16,70,60 

15 

I  r 

(HGP-7ROOO. ) 

16,80,70 

16 

lc 

(HGP-53000. I 

1  7,90,80 

17 

IF 

( HGP-47000 . ) 

18,100,90 

IB 

IF 

( HGP— 2  5000 • ) 

19,110,100 

19 

IF 

(HGP-1 1000. ) 

120,120,110 

70 

A  1 

=  .222129F-0R 

A2 

=  9.R6137 

TR 

=  2836.188 

PR 

=  2 • 97  5  9F— n6 

RHOR  =  6  # 1 1 3E-1 3 

HR 

=  200000. 

Gn 

to  no 

RO 

A 1 

=  .3507i6F-05 

A  2 

=  6.83296 

TR 

=  2566.188 

PB 

=  5.8954E-06 

RHflR  =  1 . 33BE-1 2 
HB  =  170000. 

GO  TO  130 

40  AT  =  .7S4341F-03 
A  2  =  3.41 648 
TB  =  2386.188 
PB  =  7.5578E-06 
RHOB  =  1.84SE-12 
HR  =  160000. 

GO  TO  130 

50  A 1  =  . 886289E-04 
A  2  =  1.70824 
TB  =  406.188 
PB  =  1.S562E-0A 
RHOB  =  2.232F-10 
HB  =  106000. 

GO  TO  1 30 


12 


FDL  TDR  64-112 
Part  I 

TABLE  1  (CONT’D) 

60  A 1  =  .241458F-04 

A2  =  8,5412  . . . . 

TB  =  298.188 
PB  =  2.1809F-03 
RHOB  =  4.261F-09 
HB  =  90000. 

GO  TO  130 
70  A 1  *  0. 

A3  *  .206234E-03 
TB  =  298.188 
PB  =  2.108E-07 
RHOB  =  4.1189F-08 
HB  =  79000. 

GO  TO  125 

80  A1  =  -• 1 59202F-O4 
A2  =  -7.59218 
TB  *  508.788 
PB  =  1.2181 
RHOB  =  1.39469E-06 
.  HB  =  53000. 

GO  TO  130 
90  A1  =  0. 

A3  =  . 120869E-0  5 
TB  «  508.788 
PB  =  2.5155 
RHOB  =  2 . 8804F-06 
HB  *  47000. 

GO  TO  125 

100  A 1  =  *  1 38466E-04 
A2  =  11.3883 
TB  =  389.988 
PB  =  51.979 
RHOB  =  7 • 76  5  E -0  S 
HB  =  25000. 

GO  TO  130 
110  A 1  =  0 • 

A3  =  . 1 57689E-03 
TB  =  389.988 
PB  =  472.73 
RHOB  =  7.062F-04 
HB  =  11000. 

GO  TO  125 

120  A 1  =  —.22  5569P-04 
A2  =  -5.25612 
TB  a  518.688 
PB  =  2116.71695 
RHOB  *  2.37692E-03 

HB  =  0.  —  ...  ... 

GO  TO  130 

125  P=PB*FXPF I -1 ,*A3* IHGP-HB 1) 

RHO  =  RHOB*EXPF t -1 • *A3* { HGP-HR ' 1 
GO  TO  131 

130  P  =  PR*EXPF (-1 . #A?*LOGF ( 1 .  + A  1  * { HGP-HB )  1  ) 

RHO  =  RH0B*EXPF(I-1.*(  1.+A2  }.|*LOGF(  1 .+A1* I HGP-HB 1)  J 

131  TM  =  TB#  ( 1 *  +  A 1  * { HGP-HR 1 ) 

VS  =  49. 020576*SGRTF ( TM ) 

IF  ( HGP-90000  *  1  134,134,132 

132  VS  =  0, 

VI 5  =  0. 

IF  (HGP-1 80000. I  135.135,133 
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TABLE  l  (CONT'D) 

133  IF  { HGP-1 200000  * )  136.213,213 

134  T  =  TM 

VIS  =  ( T*SQRTf(T} “.O226908E-O6 ) / I ( T+198.72 ) +RHO ] 

GO  TO  140 

135  A  =  ,739511 

_ _R  =  .174164  _ _ _ 

C  =  220. 

0  =  25. 

GO  TO  137 

136  A  =  .935797 
B  =  .273966 

C  =  J.B0. _ 

D  =  140. 

137  T  =  TM#(  A-8*ATANFM  HGP-CI  / 1>  )  ) 

140  IF  (VS)  141,141,142 

141  XMN(N)=0. 

GO  TO  143 

142  XMN  (N  )  =VA  (  N  )  /VS  . 

ALFH0=07+ALFD| N ) 

IF  (ALFHD)  143,143,144 

143  °2=P 
RH02=RH0 
T2IN)  =  T 
V2  =  VAIN) 

GO  TO  147 

144  ALFHR  =  ALFHD/57.295 

CP=2.*ALFHR**->*(  ,6+SQRTF  <  .36  +  1  ./  f  (XMN  t  N  )*ALFHR  1**21)) 

P2*P*( .7*CP*XMNf N )**2+l. I 
RH02=RF0* (6.+o2/P+l.)/(6.+P2/P) 

T?  (N ) = T *P  2*RHn / ( P*RHQ2 ) 

XMNN=SGRTF(  I6.*P2)7(7.#Pi+1«/7» ) 

IF(XMNN/XMNtN|-l . ) 1 46 *  1 46 ,  I  4 3 

146  BETA=ASIN (XMNN/XMNI N) 1 

XMN2*( l./SINF (BETA- ALFHR )  ) +SQRTFI 5./ I  6 . *RH02 /RHO- 1 . )  ) 

V2=49, ] +XHN2+SQRTF I T  ?  I  N  ) ) 

147  BMU2  =  1. 1BE-05+716.+I T?(N 1/500. )*SQRTF( T2fN) /500, )/(T2 ( Nl+216. ) 
RN2 IN ) =RH02*V?*ELH*32. 174/BMU2 


IF 

f  RN2 ( N l-RNCR)  148  »1 49, ]49 

149 

CH 

=  .3  32 

YH 

“  .5 

RH 

=  .88 

GO 

TO  1  80 

149 

CH 

=  .0296 

YH 

*  .a 

RH 

=  .9 

nn 

D1 

=  —94 .38 

r>2 

=  .2331 

D3 

-  8.4E-06 

H  2 

=  T2 l N ) **2*D3+T? { N )*D2+D1 

HAW 

“  V2  *  *?*RH/ 8  *  0 1  2  F  +  04  +H2 

IF( TS(N)-TTS I ?0)> 151,151 ,217 

151  I  *  I) 

152  I  =  1+1 

IF  (TS(N)-TTSt II 1  1  54,1  5  3,1  52 

153  ES  =  ES ( I) 

GO  TO  155 

156  FS  =  FSll-1  )  +  [TS(N)-TTS(  1-1  )  )*(ESt  I  1-FS(I-1  )  ) /(TTSt  I)-TTSU-1  )  ) 

155  J  =  0 

156  J  =  J+I 

IF  (TS(N)-TTStJI)  158,157,156 

157  CPS  =  CPSfJ) 

GO  TO  159 
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TABLE  1  (CONT’D) 


158  CP5  »CPS( J-l)+(T5(N)-TTSt J-l> ) * ( CPS f J I -CPS ( J- 1 I  > / ( ITS ( J 1 -T TS f J- 1 ) ) 

159  HS  =  T5fN }+*2*D3+T5(N(*D2+Dl 

HR  =  .22* [HAW-H2)+.5*(HS-H2 )+H2 
D4  =  400. 

D  5  =  8.829 
06  =  -1 • 978E-04 
T  A  W  (  N  )  =  HA  W**2*D6+H A W*P  5«  D4 
TH  =  HR**2*D6+HR*D5+D4 

BMUH  =  1. 18E-05*716.* ( TH/500. ]*SQRTF (TH/500. )/( TH+216. ) 

RHOH  =  P2/{ 53.3#TH) 

RNH  =  RH0H+V2  *ELH / BMUH 

SBC  =  4.758E-13  . .  .  „  _ 

DT5  (  N  )  =  ICH*RMUH+EXPF(YH*LOGF  I  RNH  )  )  *P  RP*  f  HA  W-H  S  1  )  /■(  ELH*RHO$*DELT  A 
1S*CPS  J-(SBC*E.S*(TSIN(**4-T2IN)**4)  ) / ( DELT AS*RHO$*CPS  1 

160  DTS I { 1 )  =  OTS [N ) 

DO  173  L=  1 , 3 

ISA  =  TS( Nl+OTSI (L)*DELTAH/2. 

_  HS  =.  TSA**2*D3+T5A+D2+DL  .  .  _ _  ... 

HR  =  .22* (HAW— H2  (+.5MH5-H2  (+H2 
TH  s  HR**2*D6+HR*D5+D4 

BMUH  =  1, 18E-05*716,* { TH/500. I  *  SORT F ( TH/500. ) / ( TH+2 1 6 . ) 

RHOH  =  P2/f33,3*TH) 

RNH  =  RH0H*V2*ELH/BMUH 

DTS I ( L+l )  =  (CH*BMUH*EXPF(YH*LOGF(RNH> ) *PRP* I  HA W-HS ) ) / I FLH#RH0S*DF 

1LTAS+CPS )  - ( SBC*CS+ ( T  SA**4-T2 IN ) **4)  ) / ( DFL T A S*RHOS*CPS I 
IF  ( L-2 )  170,171,172 

170  DTA  =  DTS I  ( L  +  1  ) 

GO  TO  173 

171  DTR  =  DTSI  (L+1  ) 

„  ,_DTSI(L  +  1)  =  2.*DTSI  (L  +  l)  .  . .  .  ...  _  . 

GO  TO  173 

172  DTC  =  DTSI (L+l ) 

173  CONTINUE 

0 ( N  )  =  DTSI N ) *CPS*RHOS*DFLTA5*3600. 

CPA  =  , 06 85 6* t 3. 5  + 1(5 526. /TS IN  I  ) **?*FxPF (  5  5  26 . / T5 I N )  ) / ( -1 .+EXPF { 5  5 
126,/TSIN) ) )**2 ) ) 

HTCIN)  =  CPA#CH*BMUH*EXPF ( YH+LOGF { RNH ) ) *PRP*3600 • / ELH 
200  TSIN+l)  =  TS.{*i)  +  (  DELTAH/6.  > *  I DT S ( N ) +2 . *DT A+2 . +DTB+DTC ( 
BC=ESA*SBC*3600. 

HMAX=HTC( 1 ) 

002041  =  2, < 

IFtHMAX-HTCf I))203»204»204 

203  HMAX=HTC( I ) 

204  CONTINUE 

T  A WMA  X  =  T  A W ( 1 ] 

D0206 I  =  2 , < 

IF ( TAWMAX-TAWI I  I  1205,206  ,206 

205  T  A WMA  X  =  T  A  W (  I ) 

206  CONTINUE 

IF ( TAWMAX-4000. (207,209,219 

207  IF ( TA WMA X -2000. (208,210,209 

208  IFITAWMAX-1000. (211 ,211 ,210 

2  09  BH  =  FXPF ( .25+LOGF ( 500./HMAX )  (*FXPF( . 25  +  LOGF { BC/2 . 0F-09 ) ) 

DX=500. /HMAX 
T  AWR  =  4000 . 

G0T0212 

210  BH  =  EX  PF I . 2  5  *L0GF (500./HMAX))*EXPF(.25*LnGF(BC/2.0F-09)l*.5  95 
DK=2,*500./HMAX 
TAWR=2000. 

G0T02 12 


15 


FDL  TDR  64-112 
Part  I 

TABLE  1  (CONT’D) 


211  BH=EXPF< .25*LOGF(5O0./HMAX> ) *EXPF ( .25*L0GF f BC/2 .0E-09 )1 *. 154 
0K=4. *500 • /HMAX 
TAWR  = 1 000 . 

21*  T,SMAX  =  TSm  .  ..  .  .  . . . . .  •  -  --  .  — 

D0225 I  =  2  *  X 

IF ( TSMAX-TS I  I  V )  224. 225. 225 

224  TSMAX=TS(T1 

225  CONTINUE 

IF<TSMAX-4000.)226.?28,219 
...  226  IFITSMAX-2000.  I  22.7 .229  *228  -  - - 

227  IF (TSMAX-1000. I  230. 230,229 

228  TSR=4000. 

G0T0231 

229  T SR =2 000 » 

G0T0231 

230  ,TSR=1000» _  ..  . . 

231  n0232I=l.K 

TSP ( I  I  a  T  S I  I 1*Q9.5/TSR 
T AWP (  I  I =  T AW (  I >  *99.5/TAWR 

232  HP  (  1  J  =HTCm*99.5/HMAX 
WR1TFOUTPUTTAPE3.220.FLH 

_ 220._F0RMAI...I  LHA.1X.5HELH  =.E11.4)  . . .  . 

222  WR I TFOUTPUTT APE3.223.BH 

223  FORMAT  (2X.5HRH  «.E11.4> 

WRITEOUTPUTTAPE3.233.DK 

233  FORMAT  (2X.5HDK  a.Fll.4) 

WRITE0UTPUTTAPE3.20] , ( THFTAI I } .H(!) ,VA( I ) * ALFD { I ) ,XMN( I ) . TS C I  I , 

_ lDTSfJ I.TAW1I ) *T2t IV*  RN2! 1 ) *  I =1 »K 1 

201  FORMAT  (3X.*8HTIME  =  SEC,6X,6HALT  =  FT,5X,7HVEL  =  FPF,5X,8HALFD=DEG,7X  .2H 
1 MN.9X  »4HTS=R  »6X , 7HOT  S=RPS»6X  *5HT  AW  =  R  »8X  *4HT2=R»8X»  3HRN2  7 1HA *  1  OF  12  ■ 
24/i  ix.mei2.4i ) 

WRITFOUTPUTTAPF3.202. [THFTAI  I ) »HTC ( I ) .0  t I > .TAWP(  I )  .HP |  I  )  ,  T SP [  I  I 
1 . 1  =  1 .K  ) 

20?  FORMAT! 1  HA . 2X ,8HT IMF  =  SEC • 7X . 3HHTC . lOX. 1HQ . lOX . 4HTAWP • 9X . ?HHP . 9X .3 
i HT  SP /  1  H  A  ♦  6  F 1 2 . 4  /  1 1 X  »  6  F 1  2  •  4  )  I 
GO  TO  5 


219  WRITFOUTPUTTAPE3.2Z1 

221  FORMAT  !23H  TAW  OR  TS  OUT  OF  RANGE) 

GOTO? 1 S 

213  WRITE  OUTPUT  TAPF  3.214 

214  FORMAT  (2»v.21HALTltUDF  OUT  OF  RANOFl 
GO  TO  215 

217  WRITFOUTPUTTAOF^.218 

218  FORMAT! 40H  SKIN  TEMPFR A  TURF  FXCFFOS  RANGE  OF  TABLE) 

215  CALL  FXIT 

216  FNO 
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7.  Error  printout  if  an  excess  altitude  is  inadvertently  entered  with  the  data 

8.  Error  printout  if  the  skin  temperature  exceeds  the  maximum  temperature  of 
the  interpolation  table 

9.  Error  printout  if  the  adiabatic  wall  temperature  (T  )  or  skin  temperature  (T  ) 

aw  s 

exceeds  4000°  R 

10.  Scale  measurements  of  the  heat  transfer  coefficient  (h),  the  adiabatic  wall 
temperature  (X  ),  and  the  skin  temperature  (Tg)  for  the  .heat  rate  computer 

The  following  information  is  required  as  data  for  the  program: 

1.  Number  of  time  points,  K 

2.  Flight  profile  information  (inserted  as  five  data  points  per  card) 

a.  Time  (THETA) 

b.  Altitude  (H) 

c.  Velocity  (VA) 

d.  Angle  of  attack  (ALFD) 

3.  Material  specific  heat  and  emissivity  versus  temperature  (20  points  per  curve, 
five  data  points  per  card) 

a.  Temperature  range  (TTS) 

b.  Material  specific  heat  (CPS) 

c.  Skin  emissivity  (ES) 

The  previous  data  is  required  only  once  for  each  flight  profile.  The  following  data  is 
required  for  each  point  under  consideration: 

1.  Beginning  skin  temperature  (TS(1)) 

2.  Prandtl  Number  raised  to  the  -(2/3)  power  (PRP) 

3.  Critical  Reynolds  Number  (RNCR) 

4.  Wedge  angle,  which  is  one-half  the  total  included  angle  (D7) 

5.  Distance  from  the  leading  edge  to  the  point  in  question  (ELH) 

6.  Material  density  (RHOS) 
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7.  Material  thickness  (DELTAS) 

8.  Time  increment  (DELTAH) 

9.  Average  skin  emissivity  for  expected  temperature  range  (ESA) 

The  results  of  the  computer  program  are  printed  in  tabulated  form  with  headings  at 
the  top  of  each  column.  The  following  results  are  printed  adjacent  to  the  corresponding 
time: 

1.  Time 

2.  Altitude 

3.  Velocity 

4.  Angle  of  attack 

5.  Local  Mach  Number 

6.  Skin  temperature 

7.  Derivative  of  skin  temperature  with  respect  to  time 

8.  Adiabatic  wall  temperature 

9.  Local  atmospheric  temperature  aft  of  the  shock  wave 

10.  Local  Reynolds  Number  aft  of  the  shock  wave 

11.  Heat  transfer  coefficient 

12.  Total  heat  transferred 

13.  Paper  tape  punch  values  for  the  adiabatic  wall  temperature,  skin  temperature, 
and  the  heat  transfer  coefficient 


The  tabulated  results  are  headed  by  a  printout  of  the  location  of  the  point  in  question, 
the  B-h  potentiometer  setting,  and  the  demand  gain.  All  results  are  given  to  four 
significant  digits. 

The  information  obtained  from  the  digital  computer  program  is  used  to  design  the 
laboratory  test  setup  for  simulating  the  aerodynamic  heating.  The  heat  transferred  and 
the  tempature  rise  rate  are  used  to  determine  power  requirements,  and  the  heat 
transfer  coefficient  and  the  adiabatic  wall  temperature  are  used  as  inputs  to  the  heat  rare 
computer. 
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APPENDIX  II 


SPECIAL  PURPOSE  HEAT  RATE  COMPUTER 

The  General  Electric  Heat  Rate  Computer  can  calculate  the  amount  of  energy  required 
to  simulate  aerodynamic  heating  of  forty  airframe  test  areas  and  controls  the  flow  of 
electrical  energy  used  to  simulate  these  effects.  The  computer  has  magnetic  drum  storage 
for  300  points  per  channel.  Actual  computing  is  accomplished  by  one  closed  loop  analog 
computer  which  is  shared  with  the  forty  information  channels  by  means  of  multiplexing. 

The  heat  rate  computer  is  designed  for  three  modes  of  operation:  (1)  the  power 
feedback  mode  (YEI);  (2)  the  derivative  mode  (dT  /dt);  and  (3)  the  temperature  versus 

time  mode  (T  /t).  In  the  power  feedback  mode  the  computer  operates  on  the  following 
equation: 

A.  ■  0,  [  h  <  T„  -T, )  -DT/  -YEl] 


In  the  derivative  mode  the  computer  operates  on  this  equation: 

Ae  =  G,  [  h  (Taw  -T,  >  -BTS4  -  cw  r  ] 

In  the  temperature  versus  time  mode  the  computer  does  no  calculation  but  operates 
only  as  a  controller,  comparing  the  test  temperature  with  the  desired  temperature 
programmed  on  the  computer  drum.  The  equation  for  this  mode  is  expressed  as 
follows: 


Ae 


(  T 


sp 


'st 


) 


where  T  is 
sp 


the  programmed  temperature  and  T  is  the  specimen  temperature. 

St 


The  heat  rate  computer  provides  the  error  signal  to  the  ignitron  units,  which  provide 
electrical  power  for  the  test  setup.  The  design  of  the  power  source  is  very  important 
for  good  computer  control  and  test  results.  The  power  source  must  provide  sufficient 
power,  but  not  an  excess  of  power.  An  excess  amount  of  power  will  usually  result  in 
poor  computer  control  and  a  higher  temperature  than  is  desired.  The  power  required 
is  based  on  either  the  greatest  temperature  rise  rate  required  or  the  power  required 
to  maintain  the  maximum  temperature.  Additional  power  is  required  to  allow  for  the 
efficiency  of  the  test  setup.  The  efficiencies  of  most  test  setups  in  this  laboratory 
range  between  approximately  20  and  50  percent.  Therefore,  a  power  availability  of 
about  two  to  five  times  the  calculated  required  power  will  usually  result  in  satisfactory 
computer  control. 
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SCALING 

In  the  power  feedback  (YEI)  and  derivative  (dTs/dt)  modes  the  computer  operates 
in  a  closed  loop  which  consists  of  a  demand  side  and  a  feedback  side.  The  demand  side 
is  represented  by  the  following  portion  of  the  design  equations: 

h  1T0W“TS  )  -  BTS4 

dTs 

The  feedback  side  is  represented  by  YEI  for  the  power  feedback  mode,  and  cwr  ^ — 

for  the  derivative  mode.  Both  sides  of  the  closed  loop  must  be  scaled  properly  for 
satisfactory  test  results.  A  schematic  diagram  of  the  power  feedback  and  derivative 
modes  is  shown  in  Figure  5. 

DEMAND  SCALING 

The  demand  scaling  is  identical  for  both  the  power  feedback  mode  (YEI)  and  the 
derivative  mode  (dTs/dt).  The  factors  involved  are  the  heat  transfer  coefficient  curve 
(h),  the  adiabatic  wall  temperature  curve  (T  ),  the  B-h  SCALE  FACTOR  potentiometer, 

clW 

the  skin  temperature  range  selector  switch,  the  thermocouple  selector  switch,  the  K 
SCALE  FACTOR  switch,  the  power  range  selector  switch,  and  the  mode  selector  switch. 

The  mode  selector  switch  has  three  positions:  (1)  El,  for  power  feedback  mode, 

(2)  dTs/dt,  for  derivative  mode,  and  (3)  Ts/t,  for  the  temperature  versus  time  mode. 

The  thermocouple  selector  switch  is  set  to  correspond  with  the  type  of  thermocouple 
used  on  the  test  specimen.  Positions  are  available  for  thermocouples  of:  Iron-Constantan; 
Chromel-Constantan;  Chromel-Alumel;  and  Platinum-Platinum,  10-percent  Rhodium. 

The  skin  temperature  range  selector  switch  has  positions  of  1000°,  2000°,  and  4000° R. 
This  switch  should  be  set  on  the  closest  range  that  exceeds  the  maximum  adiabatic  wall 
temperature  stored  on  the  drum  or  the  maximum  skin  temperature,  whichever  is 
greatest. 

The  K  SCALE  FACTOR  switch  has  settings  of  1,  10,  and  100  available.  Although 
this  switch  is  usually  left  at  a  setting  of  1,  It  may  be  used  to  obtain  more  demand  gain 
than  is  available  in  the  inherent  gain  obtained  in  the  scaling  of  the  h  and  T  curves. 

The  El  RANGE  switch  has  settings  that  represent  an  electrical  load  of  72,  144,  360, 
and  720  KW.  In  the  derivative  mode  (dTs/dt)  this  switch  has  no  control  feature  and  It 
should  be  set  at  the  nearest  range  that  exceeds  the  connected  load  when  the  full  voltage 
(580  volts)  is  applied.  Each  power  range  also  has  a  gain  associated  with  it.  These  gains 
are  important  in  the  power  feedback  mode  (YEI),  where  the  El  RANGE  switch  is  a  part 
of  the  feedback  mode  circuitry.  The  manipulation  of  this  switch  in  the  power  feedback 
mode  will  be  discussed  in  the  section  on  feedback  scaling. 
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DEMAND 


I - 1 


FEEDBACK 


Figure  5.  Schematic  Diagram  for  Power  Feedback  and  Derivative  Modes 
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One  heat  transfer  coefficient  curve  per  channel  is  stored  on  the  magnetic  drum  after 
being  inserted  by  paper  tape.  The  curve  is  scaled  according  to  the  following  equation: 

h  tape  value  =  99.5  ^  J 

"max 


The  IBM  digital  computer  program  performs  this  scaling  operation  with  the  exception 
of  changing  the  third  digit  of  each  point  to  “0”  or  “5”.  The  third  digit  of  each  point 
must  be  rounded  off  to  0  or  5  before  the  heat  rate  computer  will  properly  accept  the 
points. 


The  heat  rate  computer  has  provisions  for  as  many  as  ten  adiabatic  wall  temperature 
curves.  These  curves  are  also  coded  on  paper  tape  and  stored  on  the  magnetic  drum. 
The  scaling  performed  by  the  IBM  digital  computer  program  follows  this  expression: 

Taw  tape  value  =  99.5  f— — — — —  ) 

range 

Again  the  third  digit  must  be  manually  rounded  off  to  0  or  5. 


The  maximum  value  of  the  heat  transfer  coefficient  and  the  Tg  range  selection 
determine  the  inherent  demand  gain,  K^,  obtained  with  the  scaling  of  h  and  T^. 
The  equation  for  determining  the  inherent  demand  gain  is  as  follows: 


K 


D 


500 

hmax 


4000 
Ts  range 


This  calculation  is  also  provided  by  the  IBM  digital  computer  program.  Best  results 
are  obtained  if  h  does  not  exceed  500  and  the  demand  gain  does  not  exceed  10. 


The  B-h  SCALE  FACTOR  potentiometer  is  a  ten-turn  potentiometer  that  is  designed 
to  allow  for  the  radiation  losses  from  the  test  specimen.  Each  turn  of  this  potentiometer 
represents  a  value  of  1.0.  The  potentiometer  setting  is  determined  by  the  following 
expression: 


-h  =( 


B  )A  1 

f  Ts  range  \ 

f  500  V 

2  oof* 

^  4000  ' 

'  '’max  > 

'/« 


The  term  "B"  in  this  expression  is  the  product  of  the  specimen  emissivity  and  the 
Stefan-Boltzmann  constant.  The  IBM  digital  computer  program  calculates  the  B-h 
SCALE  FACTOR  potentiometer  setting  using  an  average  emissivity  over  the  expected 
temperature  range  of  the  specimen. 

FEEDBACK  SCALING 

Feedback  scaling  in  the  power  feedback  (YE*)  niode  requires  proper  manipulation 
of  the  El  RANGE  switch  and  the  K  SCALE  FACTOR  potentiometer.  Feedback  gain  is 
obtained  through  the  use  of  the  power  range  selector  switch.  The  power  ranges  of  72, 
144,  360,  and  720  KW  have  associated  gains  of  10,  5,  2,  and  1,  respectively.  Therefore, 
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the  required  gain,  as  well  as  the  required  power,  must  be  considered  in  arriving  at  the 
power  setting  of  the  El  RANGE  switch.  The  utilization  of  this  switch  will  be  discussed  in 
the  scaling  examples. 

The  K  SCALE  FACTOR  potentiometer  is  a  ten-turn  potentiometer  with  each  rotation 
of  the  dial  representing  a  value  of  0,1,  This  potentiometer  is  used  to  attenuate  the 
feedback  signal  to  allow  for  the  efficiency  of  the  test  setup.  When  the  power  range 
selector  switch  is  set  so  the  feedback  gain  is  the  same  as  the  demand  gain,  the  setting 
for  the  K  SCALE  FACTOR  potentiometer  is  the  efficiency  of  the  channel  (rj)  divided 
by  the  area  (A)  covered  by  the  channel.  However,  the  El  RANGE  switch  and  the  K  SCALE 
FACTOR  potentiometer  may  be  combined  in  such  a  way  that  almost  unlimited  values  of 
feedback  gain  may  be  obtained.  The  combination  of  these  two  controls  is  covered  in  the 
scaling  examples. 

In  order  to  arrive  at  the  proper  K  SCALE  FACTOR  potentiometer  setting,  the  efficiency 
of  the  test  setup  must  be  known.  This  efficiency  is  usually  determined  by  a  constant 
power  efficiency  test.  This  test  may  be  run  on  a  simulated  test  setup,  but  should  be 
conducted  on  the  actual  test  setup,  if  feasible.  The  efficiency  test  should  be  designed  to 
require  a  constant  power  that  will  give  a  specimen  temperature  rise  rate  that  is 
comparable  to  that  expected  in  the  actual  simulation  test. 

The  efficiency  test  is  conducted  in  conjunction  with  the  data  processing  system. 

The  data  section  has  a  program  written  for  the  CDC  1604  digital  computer  based  on  the 
following  equation: 


V 


A 


/  \  /  AT» 
'  EX  '  '  AT 


(  1.0542  ) 


This  equation  is  legitimate  as  long  as  it  is  applied  to  a  portion  of  the  temperature  curve 
where  the  rise  rate  is  relatively  constant.  To  prepare  the  efficiency  program,  the  data 
section  must  be  provided  the  following  information:  material  thickness  (t),  material 
density  (w),  a  curve  of  material  specific  heat  versus  temperature,  thermocouple  type, 
and  the  power  range  to  be  used.  The  1604  computer  is  then  operated  on-line  with  the 
heat  rate  computer  to  provide  the  temperature,  temperature  rise  rate,  power  absorbed 
by  the  specimen,  total  power  required,  and  the  efficiency  per  unit  area.  The  efficiency 
per  unit  area  is  then  used  in  setting  the  K  SCALE  FACTOR  potentiometer. 

Three  examples  of  using  the  El  RANGE  switch  and  the  K  SCALE  FACTOR  potenti¬ 
ometer  in  the  feedback  circuit  of  the  power  feedback  mode  (YEI)  will  now  be  shown. 

Only  the  feedback  scaling  will  be  given  here.  Examples  of  demand  scaling  are  shown  in 
the  discussion  of  derivative  mode  (dTs/dt)  scaling. 

Example  1: 


Power  required  =  300  KW 

7)/a  =  .312  (  dele  rmined  from  on  efficiency  test) 
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In  this  case  the  El  RANGE  switch  should  be  set  on  a  power  range  of  360  KW  since  this 
is  the  nearest  position  available  that  exceeds  the  power  requirement.  The  360-KW 
position  also  has  an  associated  gain  of  2,  which  in  this  case  matches  the  demand  gain. 
Therefore,  for  this  example,  the  K  SCALE  FACTOR  potentiometer  would  be  set  at  the 
value  of  tj/A,  or  .312.  This  value  is  3.12  turns  of  the  potentiometer  dial. 

Example  2: 

Given!  Kq  =  2 

Power  required  -  100  KW 

77/A  =  .247  (determined  from  on  efficiency  test  > 

In  this  example  the  El  RANGE  switch  may  be  placed  at  144  or  360  KW.  The  144-KW 
range  would  be  sufficient,  but  the  associated  gain  of  5  does  not  match  the  demand  gain. 

If  the  360-KW  position  were  used,  the  gains  would  match  and  the  K  SCALE  FACTOR 
potentiometer  would  be  set  at  the  value  of  v/A.  Optimum  utilization  of  the  full  power 
range  would  be  obtained  with  the  144-KW  setting.  However,  the  K  SCALE  FACTOR 
potentiometer  setting  must  be  reduced  to  compensate  for  the  difference  between  the 
demand  gain  and  the  feedback  gain.  This  can  be  done  by  multiplying  the  17/ A  value  for 
matched  gains  by  the  ratio  of  the  demand  gain  to  the  gain  associated  with  the  chosen 
power  range.  The  new  setting  for  the  K  SCALE  FACTOR  potentiometer  in  this  case  would 
then  become: 


Kpot  "  V/A  (  EI  9<jin  ) 

■  <-2471  (t) 

K  p0t  =  .0988 


Example  3: 

Given  !  K q  -  4.  72 

Power  required  ~  100  KW 

^/A  =  .417  (determined  from  on  efficiency  test  ) 

In  this  example  the  gain  value  is  such  that  no  power  range  selection  will  match  it.  A 
power  range  of  144  KW,  with  an  associated  gain  of  5,  may  be  used.  The  rj/A  value  for 
matched  gains  must  then  be  modified: 

Kpot  '  41 '  (ex  goin  / 


Kpot  a  -394 
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A  general  expression  for  finding  the  proper  K  SCALE  FACTOR  potentiometer  setting 
can  be  derived  from  the  previous  examples.  This  expression  is: 


where  El  gain  is  the  associated  gain  of  the  chosen  power  range.  This  equation  may  be 
used  as  long  as  the  calculated  value  does  not  exceed  1,0.  For  instance  if  the  demand 
gain  is  5  and  the  tj/A  is  .25,  the  El  RANGE  switch  may  be  set  at  360  KW,  or  a  gain 
of  2.  The  potentiometer  setting  would  then  become: 


Kpot  =  t  •  25  )  (  2  ) 
Kpof  =  .  52 5 


However,  under  these  conditions  the  El  RANGE  switch  could  not  be  placed  at  720  KW,  or 
a  gain  of  1,  because  the  potentiometer  setting  would  then  become: 


Kpot  s  <  (t) 


Kpot  =  1.25 


which  is  impossible  to  set  on  the  potentiometer.  Therefore,  the  power  range  and  the 
demand  gain  must  be  combined  in  a  way  such  that  the  K  SCALE  FACTOR  potentiometer 
calculation  does  not  exceed  1.0. 

The  feedback  scaling  for  the  derivative  mode  (dTs/dt)  is  very  similar  to  that  of  the 
power  feedback  mode.  In  the  derivative  mode  the  C  SCALE  FACTOR  switch  and  the  C 
SCALE  FACTOR  potentiometer  are  used  to  obtain  the  required  feedback  gain.  The  C 
SCALE  FACTOR  switch  is  a  four-position  switch  with  gain  settings  of  1,  2,  5,  and 
34.125*.  The  C  SCALE  FACTOR  potentiometer  is  a  ten-turn  potentiometer  with  each 
rotation  of  the  dial  representing  a  value  of  1.0.  When  the  C  SCALE  FACTOR  switch  is 
set  on  a  gain  that  matches  the  demand  gain,  the  proper  setting  for  the  C  SCALE  FACTOR 
potentiometer  is  cwr.  However,  the  C  SCALE  FACTOR  switch  and  the  C  SCALE 
FACTOR  potentiometer  can  be  combined  to  give  almost  any  value  of  feedback  gain.  The 
use  of  these  two  controls  is  discussed  in  the  scaling  examples. 


♦This  gain  is  obtained  by  setting  the  C  SCALE  FACTOR  switch  at  10.  The  switch  was 
originally  designed  with  a  gain  of  10,  but  a  resistor  was  changed  during  computer 
modifications  to  give  the  new  gain  of  34.125.  This  position  is  now  intended  for  the 
temperature  versus  time  mode  only,  but  it  may  still  be  used  in  the  derivative  mode. 
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Two  examples  of  derivative  mode  scaling  will  now  be  given. 


Example  4: 

Given!  hmox=  200  BTU/ hr -ft2-  *R 
Ts  range  =  2000“  R 
K  SCALE  FACTOR  switch  =  I 
€  =  *6 

cr  =  1.73  X  I0-9  BTU/hr-ft2-  *R4 
c  =  .  12  8 T U  /lb  -  *R 
w  =  500  lb/ ft3 

t  =  .06  in.  =  ,005  ft 

First  determine  the  inherent  demand  gain,  K^: 

,  500  v  ,  4000  ^ 

D  '  hmo*  Ts  range 

/  500  w  4000  \ 

V  200  / V  2000  / 


Now  determine  the  setting  for  the  B-h  SCALE  FACTOR  potentiometer: 


f  B  f* 

/  Ts  ronge 

)A 

/  500 

'  2  X  10"*' 

'  4000 

) 

'  hma> 

€  cr 

(1.73  X  10"*  ) 

(.6) 

1.038  X  10"* 

/  1.038  X  lO-*  ■ 

k'4  /  2000 

f4  ( 

500  y 

'  2  x  10"* 

/  l  4000 

}  \ 

200  ' 

=  (.  si  s)/4  ( .  500 ) iz.5)^ 

=  (.849  J  (.595)  (1.258  ) 

B - h  =  . 635 


Now  the  feedback  side  must  be  scaled.  Since  the  demand  gain,  K^F  is  5,  the  C  SCALE 

FACTOR  switch  can  conveniently  be  set  on  a  value  of  5.  Under  these  conditions  the  C 
SCALE  FACTOR  potentiometer  is  set  at  a  value  of  cwr. 
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cw  r 

(12  1  (  500  )  (.005  ) 

.300 

Example  5: 

Given:  hmo*  =  240  BTU /hr  —  ft2  - *R 

Tsronge  =2°00°R 
K  SCALE  FACTOR  switch  =  I 
€  =  .6 

or  =  1.73X10" 9  BTU/hr  — ftZ-  •R4 
c  =  .12  BTU  /  lb  -  *R 

w  =  500  lb /ft 5 
r  =  .06  in.  =  .005  ft 

The  calculations  for  this  example  are  as  follows: 

/  4000  \  /  500  \ 

D~  \  2000  )  \  240  ) 

K0=  4.17 

B  =  ear  =  (.6)  (  1.73  X  I0~*  ) 

B  =  1.038  X  I0~a 


(  B  V/4 

f  Tt  range  , 

500  , 

f/4 

\  2  X 10”*  / 

l  4000  /  l 

hmat 

) 

j  1.038  XiO"* 

\'/4  f  2000  \% 

/  500 

-V4 

2  XIO"* 

'  \  4000  / 

v  Z40 

) 

=  (.5I9)'/4  t.500)3^4  (2.08)^ 

=  (.849)  (.595  )  (  1.201  ) 

B-h  =  .607 

cwr  =  <  .12  )  (  5  00  )  (.005) 
cw  t  =  , 300 

In  this  case  there  is  not  a  convenient  setting  of  the  C  SCALE  FACTOR  switch  that  will 
match  the  demand  gain  of  4.17.  However,  the  C  SCALE  FACTOR  switch  may  be  set  at  5. 
The  cwr  calculation  must  now  be  modified  so  the  feedback  gain  will  match  the  demand 


cpo  t  = 
cpot 
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gain.  This  can  be  done  by  multiplying  cwt  by  the  ratio  of  the  demand  gain  to  the  C  SCALE 
FACTOR  switch  position.  Therefore,  the  setting  for  the  C  SCALE  FACTOR  potentiometer 
will  be: 


'pot 


Cpot 


:w  r  ( 


kd 


Cswi tch 

«■«  (- t I) 

.250 


The  C  SCALE  FACTOR  switch  could  also  be  set  at  1.  In  this  case  the  C  SCALE 
FACTOR  potentiometer  setting  would  be: 

Cpo,  =  (.30  )  (  ) 

Cpot  «  1.231 


From  the  two  previous  examples  the  following  general  expression  can  be  derived 
for  determining  the  setting  for  the  C  SCALE  FACTOR  potentiometer: 


Cpot 


CWT 


kd  x 

^switch  ' 


This  expression  may  be  used  as  long  as  the  calculated  setting  does  not  exceed  10. 
TEMPERATURE  VS.  TIME  MODE 


In  the  temperature  versus  time  mode  the  heat  rate  computer  is  used  only  as  a  controller. 
The  error  signal  is  obtained  from  the  difference  between  the  programmed  temperature  and 
the  specimen  temperature.  The  error  signal  is  proportional  to  the  amount  the  specimen 
temperature  drops  below  the  desired  temperature.  The  following  proportional  bands  are 
available  for  the  appropriate  temperature  range: 

Scale  Proportional  Bands 

1000°R  10,  20,  30,  40,  50,  60,  70,  80,  90,  100°R 

2000°R  20,  40,  60,  80,  100,  120,  1.40,  160,  180,  200^1 

4000°R  40,  80,  120,  160,  200,  240,  280,  320,  360,  400°R 

Each  proportional  band  figure  means  that  a  specimen  temperature  lag  of  that  value  will 
result  in  the  full  error  signal  of  50  volts. 

The  power  arrangement  is  extremely  important  in  the  temperature  versus  time  mode. 
The  temperature  program  should  be  designed  so  that  no  step  increment  will  require 
more  than  one-fourth  of  the  total  available  load.  This  means  that  for  smooth  computer 
control  the  ratio  of  the  power  output  per  volt  of  Ae  should  be  held  to  a  minimum.  A 
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proportional  band  should  be  chosen  such  that  the  maximum  step  increment  will  demand 
not  more  than  one-half  kilowatt  per  lamp.  If  this  proportional  band  is  unsatisfactory, 
some  trial  runs  may  be  necessary  to  decide  which  proportional  band  gives  the  best 
combination  of  good  test  results  and  smooth  control. 

The  demand  signal  in  the  temperature  versus  time  mode  is  obtained  in  the  scaling  of 
the  specimen  temperature  by  the  digital  computer  program.  The  skin  temperature  is 
scaled  according  to  the  following  equation: 

Tt  top*  volu«  =  99.5  (  % - ) 

As  in  the  previous  cases  the  third  digit  of  each  point  must  be  rounded  off  to  0  or  5. 

The  only  controls  in  the  feedback  circuit  of  the  temperature  versus  time  mode  are 
the  C  SCALE  FACTOR  potentiometer  and  the  C  SCALE  FACTOR  switch.  With  the  K 
SCALE  FACTOR  switch  in  the  demand  side  set  at  10,  and  the  C  SCALE  FACTOR  switch 
set  at  a  gain  of  34.125*,  the  computer  technicians  set  the  C  SCALE  FACTOR  potentiome¬ 
ter  so  that  a  temperature  lag  equal  to  the  chosen  proportional  band  will  result  in  an  error 
signal  of  50  volts.  This  is  accomplished  by  first  inserting  the  maximum  desired  tempera¬ 
ture  in  the  demand  circuit.  A  millivolt  signal  is  then  manually  inserted  in  the  feedback 
circuit  to  represent  the  thermocouple  signal  equal  to  the  maximum  desired  temperature 
minus  the  chosen  proportional  band.  The  C  SCALE  FACTOR  potentiometer  is  then 
adjusted  until  the  error  signal  is  50  volts.  The  demand  and  feedback  circuits  are  now 
considered  ‘‘balanced.'’ 

To  eliminate  any  lag  in  the  temperature  versus  time  mode,  the  amount  of  demand  must 
be  increased.  This  may  be  done  in  one  of  two  ways.  (1)  each  point  of  the  demand  curve 
may  be  multiplied  by  this  factor: 

T,  m0,  +  P-  B.  Off  sat 
T*  max 

Or,  (2)  the  C  SCALE  FACTOR  potentiometer  setting  at  balanced  conditions  may  be 
divided  by  this  factor: 


_ T»  mox _ 

T.  mQx  *“  P-  B.  Offset 

The  proportional  band  offset  may  be  determined  from  the  power  arrangement.  However, 
if  the  calculated  value  does  not  give  satisfactory  test  results,  a  new  proportional  band 
offset  can  be  easily  inserted  by  using  the  second  method. 


*This  gain  is  obtained  by  setting  the  C  SCALE  FACTOR  switch  at  the  indicated  value  10. 
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The  use  of  these  two  factors  will  be  shown  in  the  following  example  of  temperature 
versus  time  scaling.  In  this  case  it  is  assumed  that  the  proportional  band  offset  is 
equal  to  the  chosen  proportional  band. 

Example  6: 

Given  '  Desired  maximum  T*  ~  I600*R 
Ts  range  switch  =  2000*  R 
Proportionol  bond  -  P.  B.  Offset  =  80®  R 

When  using  the  first  method,  each  point  of  the  temperature  curve  would  be  multiplied 
by  the  following  factor: 

T»  max  +  P-B.  Offset  _  )6QO  4  80 
Ts  max  1600 

1680 

1600 

=  I.  05 

In  preparing  the  computer  for  this  method,  the  technician  would  insert  a  voltage  in  the 
demand  side  that  would  represent  1680°R,  and  a  voltage  representing  1600°R  would  be 
applied  to  the  feedback  circuit.  The  C  SCALE  FACTOR  potentiometer  would  then  be 
adjusted  so  the  error  signal  would  be  50  volte. 

Suppose  for  illustration  purposes  the  C  SCALE  FACTOR  potentiometer  setting  was 
.495.  If  the  second  method  was  used,  the  C  SCALE  FACTOR  potentiometer  setting  would 
be  divided  by  this  factor: 

T,  max  _  1600 

Tsmax  -  R  B.  Offset  ”  1600-  80 

1600 

1520 

=  1.051 


The  potentiometer  setting  would  then  become: 

.495 

CPot  "  I  .  053 
0  pot  =  470 

By  using  the  second  method,  the  original  temperature  curve  as  Bcaled  by  the  digital 
computer  program  can  be  used  without  alteration.  This  method  also  provides  a  method 
of  easily  changing  from  one  proportional  band  offset  to  another  if  a  trial  and  error  system 
is  used  to  determine  the  best  proportional  band  offset. 
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